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Abstract. 

^ I In this article, I have investigated statistical mechanics of a non-stretched elastica in two 

dimensional space using path integral method. In the calculation, the MKdV hierarchy naturally 
appeared as the equations including the temperature fluctuation. I have classified the moduli of 
the closed elastica in heat bath and summed the Boltzmann weight with the thermalfluctuation 
' over the moduli. Due to the bilinearity of the energy functional, I have obtained its exact partition 

function. By investigation of the system, I conjectured that an expectation value at a critical point 
I of this system obeys the Painleve equation of the first kind and its related equations extended by 

■ the KdV hierarchy. Furthermore I also commented onthe relation between the MKdV hierarchy 

i and BRS transformationin this system. 

§1. Introduction 

Elastica has sometimes appeared in the history of mathematical physics according to refs.[l-4]. 
Q I The problem of elastica, an ideal thin elastic rod, was proposed by James Bernoulli in 1691. By 

c/3 ■ investigation on the behavior of an elastica, Bernoulli's family and their related people, Euler, 

dAlembert and so on, discovered many non-trivial mathematical and physical facts, e.g., classical 
field theory, minimal principle, elliptic function, mode analysis, non-linear differential equation 
and others [1-4]. In fact, James Bernoulli derived the elliptic integral related to the lemniscate 
function in 1694, before Fagnano found his lemniscate function [1-3], and found that the force of 
elastica is proportional to inverse of its curvature radius [1]. His nephew, Daniel Bernoulli followed 
James's discoveries and discovered the energy functional of elastica and its minimal principle 
around 1738. Succeeding Daniel's and James's discoveries, Euler derived elliptic integral of general 
modulus as a shape of elastica using Daniel's minimal principle and numerically integrated it. 
Then he completely classified shapes of static elasticas by numerical sketch [1], which are, nowaday, 
known as special cases of loop soliton [5]. In the computations, Euler used the static sine-Gordon 
equation. These computations essentially imply discovery of the integrable nonlinear differential 
equation, the investigation of its moduli and first application of algebro-geometrical functions to 
physics. It should be noted that they came from the studies on the elastica. Furthermore it is 
well-known that the elastica problem is the simplest prototype cj-model — > S0(2) or S0(3) 
[6,7], which was found in 18th century and investigated by Kirchhoff in last century [4]. 

Thus the elastica problem sounds to be legacy before last century but I believe that its prop- 
erties are not completely understood and its role in mathematical physics is still important. The 
difficulty to solve the elastica problem is that one must consider the constraint condition such 
as the isometry (non-stretch) condition and boundary condition. In fact, even though it is not 
an elastica problem, Goldstein and Petrich naturally rediscovered the MKdV hierarchy through 
(virtual) dynamics of a space curve with isometry condition [8,9]. (Readers should noted that the 
time-development of the physical elastic rod is not governed by the MKdV hierarchy in general 
[7,9-14] even though in some papers it seems to be misunderstood.) Due to their work, the MKdV 
hierarchy can be geometrically realized. Furthermore, using their construction, I proved that the 
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geometrical language of a space curve [11]. The classical analogue of the vertex operator naturally 
appears as a complex tangent vector of the elastica [11]. 

However these approaches are just mathematical and geometrical ones but are not directly 
connected with physical problem because in most of works there is little physical reason why the 
MKdV hierarchy must appear in physical system [8,9]. 

On the other hand, the study of elastic chain model of deoxyribonucleic acid (DNA) is current 
[16-22]. DNA usually occurs as a double helix with two complementary nucleotide chains winding 
around a common axis. In cells, the common axis of the looped DNA often folds into intricate 
structure, or a supercoiled form [16,17]. Due to the enormous size of DNA, conformation of double 
helical polymer needs topological and geometrical studies [23]. Thus the double helical polymer 
is often modeled as a thin elastic isotropic rod or an elastica and studied from the kinematic and 
topological point of view. Topological invariance such as the linking number classifies the shape 
of the DNA in three-dimensional euclidean space M"^ [16]. Based on the Kirchhoff's model of 
an elastica in and the nonlinear Schrodinger equation, the possible kinematic conformations 
of DNA were considered [16-21]. Partial thermal effect on their shapes were argued in ref.[18]. 
Furthermore the molecular dynamics study combining with the elastic energy model was reported 
[22] and then an topology change related to the knot was found. 

Statistical mechanics of elastic chain as a model of a large polymer like DNA was studied by 
Saito et al. using the path integral [24]. They calculated some exact partition function of the 
energy functional of elastic chain. Such chain model is sometimes called wormlike chain model [25]. 
Their approach is successive in the polymer physics and influences recent works [26 and references 
therein]. However they did not pay any attention upon isometry condition in a calculation of 
the path integration even though they imposed isometry condition after computation (a kind of 
quantization); they summed allover configuration space without isometry condition. It should be 
note that the constraint does not commute with quantization and statistical treatment [27]. (An 
example of such inconsistency is shown in ref.[28].) 

It is a natural assumption that flexible polymers which cannot stretch as a classical object (at 
zero temperature) cannot stretch even in heat bath. Hence it is very important to calculate the 
partition function of the elastica with isometry condition and sum the Boltzmann weight over 
only allowed conformations. 

The purpose of this article is to investigate the behavior of a closed large polymer like DNA 
in heat bath of two-dimensional space, whose shape is determined by its elasticity and stretch is 
negligible, and to clarify the physical meaning of the MKdV hierarchy. 

Furthermore it is known that DNA sometimes exhibits topology changes and has topological 
isomer [22,29]. In this article, as I will deal with an elastica in two-dimensional space, there is 
no knot invariance but exists writhing number as a topological invariant of the elastica in two- 
dimensional space if crossings are allowed [28]. Thus in this article, even though the elastica in 
two-dimensional space will be investigated, I will allow the crossings if it can be realized when I 
embed it in three-dimensional space. Then I will argue the topological change on this problem. 

The organization of this article is as follows. §2 reviews classical shape of an elastica in two- 
dimensional fiat space adding infinity point, i.e., CP^. In §3, I will investigate the statistical 
mechanics of an elastica. First, I will consider the thermal fluctuation of the extremum point 
of a partition function of the elastica in terms of the path integral method. Then I will obtain 
the MKdV hierarchy by physical requirements. Second, I will investigate the moduli space of the 
quasi-classical elastica, which is the extremum point of the partition function. Finally I sum the 
Boltzmann weight over the moduli space and obtain an exact formulation of the partition function. 
In §4, I will discuss the obtained partition function and I will comment upon the relation between 
the Goldstein-Petrich method and the BRS quantization of the gauge field [31] and a critical point 
of this model. 
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Here I will quickly review a shape of a closed elastica in two-dimensional space [10-12,32]. I 
will denote by C a shape of the elastica embedded in the projective complex line (or the Riemann 
sphere) CP^ and by X{s) its afhne vector : 

3 X{s) eC CCP\ X{s + L) = X{s), (2-1) 

where L is the length of the elastica. The Prenet-Serret relation will be expressed as 

exp(i(/)) = dsX, 9<j exp(i(/)) = i/cexp(i(^), (2-2) 

where (p is a real valued function of s and k is the curvature of the curve C, k := ds<p'- 4>{s + L) = 
(p{s) and k{s + L) = k{s). Here I have chosen the metric of the curve as the induced metric from 
the natural metric of C C CP^; by the choice (f) is real valued. 

As Daniel Bernoulli suggested to Euler [1], the energy integral of the elastica is given as 

E = [ dsk'^, (2-3) 

and shape of the elastica is realized as its stationary point. Here I assume that the elastica does 
not stretch and preserves its infinitesimal length; since deformation of the elastica is regarded 
as one parameter transformation, the assumption implies that the transformation is isometric 
[10-13,32]. Thus I will refer this condition as isometry condition. 

Following the minimal principle, I will derive the differential equation. I will consider the 
variation of the curve C — > Ce under isometry condition [8-11], 

X^X, = X + £(C/i + i[/2)exp(i<^), C/i(L) = C/i(0), C/2(L) = C/2(0) (2-4) 

where eUi and eU2 are infinitesimal real valued functions. 
Since the infinitesimal length of the curve is given as 

ds^ = dXdX = d.Xd.Xds'^, (2-5) 

corresponding length of the C^, 

dXedXe = (1 + e (ids - ik)Ui - \{ds - ik)U2)) (1 + e {{d, + ik)Ui + i{d, + ik)U2)) ds^ 

= (1 + 2e {dsUi - kU2)) ds^ + 0(e2), (2-6) 

must be ds^ modulo due to the isometry condition. Hence I obtain the relation, 

dsU^ = kU2. (2-7) 

The tangential angle of Ce is given as 

= TlogS^Xe 
1 

= (l) + e{k + dsk-^ds)Ui. (2-8) 

Its curvature is expressed as 



ke = exp(-i<^,)afX, 
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Finally I obtain the variation of the energy functional 

j\skl = j\s (^^ + £ {^d^ik") + d,^^ +0{e^). (2-10) 
Prom the variational equation, 

HI-. 

the non-linear differential equation is given as the equation of the shape of the static elastica, 

5.(^A:2 + ^) = 0, (2-12-a) 

and thus 

aifc+ + = 0, (2-12-b) 

where oi is the integral constant. This equation is known as the static MKdV equation in the 
soliton theory if derivate it by s again. 

First I will note that (2-12)'s are also equations of the energy functional [10,11], 

E= f ds (A:^ +^icos(/) + A2) , (2-13) 
Jo 

where the second term means the constraint for the relative position of X(0) — ^{L) and the 
third one is for the total length L [7,10,11,13]. The sufficiency of the third term is trivial. From 
(2-9), the second term becomes 

Ai j dscos(0e) = Ai J ds (^cos(0) - e (^kU^ + ds^^^ sin((/))j + 0{£^) (2-14) 

and by the partial integration, the second term in rhs vanishes for any Ui. This vanishing occurs 
owing to the compatibility between the MKdV equation and the static sine-Gordon equation, 

af0 + Aisin(/) = O. (2-15) 

(2-15) comes from natural variation of (p of (2-13) [1,4,7,10,11,13]. Hence (2-12)'s can be also 
regarded as the stationary point of (2-13). In fact it is known that solutions of (2-12) are in 
agreement with those of (2-15) as I will show later. It should also be noted that (2-13) is the 
(7-model with the topological term and was discovered in 18th century. In other words, the system 
of the elastica can be regarded as S0(2)-principal bundle over and the cosine term in (2-13) 
is a local version of the fundamental group vri(S'^) = Z [6]. 

Solutions of (2-12) are completely expressed by the elliptic functions. Multiplying dgk and 
integrating s [4,11,13,14,32,33], I obtain the relation, 

{d,kf = ~{k^ - ai^ + 02). (2-16) 



Introducing the quantities, P2 — P2 '■= oi, I = Pi/{Pi + P2), and P3 := -v/^r+~|32/2, k is 
expressed by the Jacobi elliptic function [34], 
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Transformation from the solutions of (2-12) to those of (2-15) is given by the identities of the 
integrand in the eUiptic integral between trigonometric function and polynomial expressions [34] . 
After integrating the differential equations, I obtain [1,4,10,11,13,33], 

2 21 
X{s) = — £;(am(/33s),0-s-i-(cn(/33s)-l), (2-18) 

P3 Pa 

where E{-, I) is the incomplete elliptic integral of the second kind and am is the Jacobi amplitude 
function [34]. Due to the closed condition, 

X{0) = X{L), (2-19) 

there is an eight-figure shape [1,4,11,13,16,33]; the modulus of the elliptic function is I = 0.908909 • • 
and the ratio of the fundamental parameters is K' /K = 0.70946 • • • . 

Thus in the set of solutions of (2.12), there are only two closed elasticas in C up to translation of 
their centroid, global rotation and scaling; circle k = 2tt/L and eight-figure shape. Here though I 
have chosen solutions such ask = 27r/L, there also exist other solutions like {k = 27rn/L | n > 1}. 

On the other hand, by taking limit L ^ oo and by considering ones in CP^, more various 
closed elasticas in CP^ are allowed. These solutions were classified by Euler in 18th century [1,4]. 

In this article the set of elasticas obtained as solutions of (2-12) will be denoted as ©cis 

©cis = {C\C is a solution of (2-12)'s}, (2-20) 
and the energy functional is expressed by, 

^cis[C] = / dsk'', Ce 6cis. (2-21) 
Jo 



§3. Statistical Mechanics of Elastica 

In this section, I will consider statistical mechanics of a closed elastica and investigate its 
behavior at heat bath. I continue to allow the crossing of the elastica even in two-dimensional 
space. I set up that there are many independent laboratory dishes in which large polymers like 
DNA individually exist one by one. A looped elastica is in the liquid whose temperature is 
determined and viscosity is very large. The liquid is a kind of heat bath. Then the kinetic energy 
of the elastica is suppressed in equilibrium state due to dissipation but owing to the fluctuation by 
temperature noise, the elastica sometimes jumps from a quasi-stable state to other quasi-stable 
states. 

The partition function of the elastica is given as [24] , 



/ 




DXexp l-P ds [/c^l . (3-1) 



Here I go on to prohibit that the local length of the elastica does change; the isometry condition 
will be maintained. In [24], it was written that if one also dealt with the kinematic term, it would 
be decoupled with the potential term (2-3). However as I employ the isometry condition, this 
statement cannot be guaranteed at all because the kinetic term is also restricted by the isometry 
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Quasi-Classical Motion. 

By the quasi-classical method in path integration [35,36], I will evaluate the partition function 
(3-1). I will expand the affine vector around an extremum point of the integral, 

X = Xqei + e(ui(s) + m2{s)) cxp(i(/.qei) + 0{e^), (3-2) 

where e is an infinitesimal parameter, e oc 1/ y/j3 and (pqci is the tangent angle of the quasi-classical 
curve of elastica. ^qci is an affine vector of the extremum point of the functional integral and will 
be determined later. In the path integral, terms with higher orders of e also play important role 
and thus I must pay attentions upon the higher perturbations of e here. Hence I will assume that 
X is parametrized by a parameter t and the difference between X and Xqd can be expressed by, 

X{s, t) := e^^*Xqei(s, t), ed^X^a = X- Xq^ + 0{e^). (3-3) 

Since for an analytic function f{x), e"^"=/(x) = f{x + a), X{s,t) can be expressed as X{s,t) = 
Xqci{s,t + e), where direction t differs from that of s; the domain of functional integration (3-1) 
deviates from the domain of the classical map (2-1). Then (3-2) becomes 

^t^qci = (""1 + i'"2) exp(i^qci), ui{L) = ui{0), U2{L) = n2(0). (3-4) 

This is virtual dynamics of the curve [8-10]. As well as the argument in §1, due to the isometry 
condition, I require [dt, dg] = for X. Then the isometry condition exactly preserves, ds = dsqd 
by the definition (3-3). This isometry relation should be compared with (2-6) which is isometry 
modulo e^. It also should be noted that although e is constant, dependence of the variation upon 
the position s is given though (3-4) and tta(s), a = 1,2. Thus (3-3) is not trivial deformation in 
general. 

I have the relation, 

-dt exp(i^qci) = {{uis - U2kqc\) + K'^2s + uik)) exp(i^qci). (3-5) 

Noting that (p ^-iid u's are real valued, 1 obtain (2-7) again from the first term and solve the 
differential equation between ui and U2 [8,10], 

dgUi = kc^clU2, = y" U2kqci =■ d~^U2kqci- (3-6) 

Here has a parameter c G M as an integral constant and coincides with the pseudo-differential 
operator. 

Then (3-5) is reduced to the equations, 

^t^qcl = -^U2, := 5f + S^feqclS^^feqcl- (3-7) 

Prom (3-3) and [dt,ds] = for X, cf) is calculated as, 

4>{S, t) = <^qcl(s, t + e)= e'^*^qcl(s, t) = (^qcl + e^tt^qd + ^6^3^ ^qd H . (3-8) 

Then noting k'^{s,t) = k^^i{s,t + e), the energy functional is expressed as 

J k^ds = J {kl^i + 2ekqcidtkqci + e^{{dtkq^cif + kc^cl^fkqcl)-\ ) ds 

= J (^qci + 2eA;qcif^ti2 + e^{(dtkqcif + ^qci^t ^qci) H )ds 



7 



Using (3-6), if I will perform the functional derivative oi E in ui, 1 obtain the classical equations 
(2-12) 's again. Since the quasi-classical configuration is realized as the extremum point of the 
functional space, I must impose the relation, 

6^^^E = 0. (3-10) 

Noting the relation (3-6), if flu2 could be regarded as another function U2 of the variation of 
the normal direction in (3-2), I might find the relation 



J dskqci^y'2 ~ J dskqciU2 = J dsdgu'i = 0. (3-11) 



Accordingly I will introduce the sequence for mathematical times t := (ti, ta, ts,- • • , t2n+i, • • ■ ) 
so that (3-11) is satisfied. I will redefine the fluctuation (3-2) and introduce infinite parameters 
family, which can sometimes become finite set as I will show later, 

X = e(Vx/^)E„=o<5*-+.a.,„^,^^^^ ^ ^^^^ ^ ^^/^^ ^ (^Wi^Wi^qci + 0(1//?). (3-12) 

n=0 

Here e was replaced with {1 / y/P)dt2n+i and 5t2„_,.i^qci is expressed as 

9t..+ ,^qcl = iui''^ + exp(0qcl), ^^"^ = 57'fcqcl4"\ 4"^ = (3-13) 

with integral constants c's vanishing for n > 1. iij'^^ is an appropriate function of s. 
Then the variation of the energy functional is calculated as, 

eds = [ {k%, + (2/v^) ^ <5t2n+lfeqcl5Wlfcqcl)ds + 0(l//3) 

n 

{kl,, + (2/V^) St2n+lKcl^ U^^^)ds + 0(1/ P) 
n 

2n+l / dsd,v!c^^^ + o{i/p) 

dskl, + 0{l/P). (3-14) 



Thus for the variations along the directions, the energy of the system is invariant modulo 
Without work, we can move it for these directions 6t2n+id2n+iX [10,37,38]. 
However for this sequence, the infinite differential equations appear [8,9], 

5wi^qci = -^^"4°^. (3-15) 

The recursion equations (3-13) are determined by the initial condition '^2°^. By following the 
thought of the quasi-classical method of the path integral, this sequence must contain the classical 
equations (2-12) 's. On the other hand, for a close elastica, there is a trivial continuous symmetry 
which is the translation of the coordinate system s along the curve C. Hence (3-15) should also 
include such translation symmetry. Of course, they contain other equations as quasi-stable shapes 
as I show following. 

Though there might be other choices, I will select minimal subspace of the functional space in 
order to satisfy above requirement. As I performed the variational computation in §1, I choose 
initial state as 
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Then the minimal set of the virtual equations of motion, which satisfies the physical requirements, 
is 

dt^„+^k^ci = -^"'d^kqcu 5*2^+3 ^ci = ^^5t2„+ifcqci- (3-17) 
First several equations and n's are given as follows, 

= 1, 4°^ = 0, (3-18-a) 

4'^=^fcqcb 4'^=5,fcqcl, (3-18-b) 

(2) 3 1 ^ ^ 3 

1 = ■^^qcl ~ 2^^*^^"^'-' ~^kdgkqc\, U2 = 2^qcl^s^qcl + t^g^qcl) (3-18-c) 



U 



77, = : dt^k^ci + dskqd = 0, (3-19-a) 

n = l: 5t3/Cqcl + 5f/Cqcl + ^fc2^l5«A:qcl = 0, (3-19-b) 

15 5 
n = 2 : St^fcqci + Sf^qci + —k^^^dsk^d + -^{dsk^^if 

+ ^fcqcl^'^cl + lOfcqcia.feqclO^feqei = 0. (3-19-c) 

Since n is identified with the Gcl'fand-Dikii operator for the MKdV equation, (3-17) is regarded 
as the MKdV hierarchy and n^^-^'s arc hamiltonian of the MKdV hierarchy [8-10,39]. 

Next I will consider the solutions of these equations. (3-19-a) means the freedom of choice of 
the origin of s and has only trivial solution k{s — ti); ti is the origin of s. Hence I must pay my 
attention only upon properties of (3-17) for n > 1. Derivative of (2-12-b) can be described as 

cd,k = ndsk. (3-20) 

By ts = s/c, this is identified with n = 1 equation (3-19-b). Here dgk of the solutions of (2-12-b) 
is interpreted as the eigenvectors of CI and c is an eigenvalue of the operator J7. Thus (3-17) 
becomes 

dt,„^,k = -n^'dsk = -c^dsk. (3-21) 

Hence the solutions of (2-12) can be solutions of all equations in (3-17). In fact for a stable 
solutions k =constant, all Uj"'' = ^'^d thus it satisfies all equations in (3-17). 

It should be noted that (3-21) comes from the fact that u^^^'^s agree with the hamiltonians 
of the MKdV hierarchy and are regarded as conservation quantities for n-th equations of n < m 
[39]. Hence using soliton theory, any solutions of (3-19-b) are solutions of higher equations in 
(3-17) n > 2. Due to the requirement of the quasi-classical solutions, any solutions must satisfy 
all of n-th equations (n > 1) in (3-17). Hence I should deal with only the solutions of the MKdV 
equation (3-19-b) as the quasi-classical solutions of this system. 

Then the sequence (3-17) fulfills the physical requirements. In other words, the solution space of 
the MKdV equation (3-19-b) is the minimal space containing the classical solutions and translation 
symmetry and filling the quasi-stable condition (3-10). 

Since for the variations along the directions t of (3-17), the energy of the system (3-14) is invari- 
ant, the fluctuation of the quasi-classical shape /cqci(s, i2n+i) should be regarded as (generalized) 
jacobi-fields or the Goldstone mode [37,38] even though they do not obey a linear differential 
equation in general. 

Here I will remark that the MKdV hierarchy naturally appears by physical requirement. It is 
very surprising because the MKdV hierarchy has infinite hamiltonians and time axes; in classical 
theory, these quantities cannot be physically interpreted. Furthermore it should be contrasted 
with the works related to space curves in [8,9,33], in which the authors chose (3-16) by hand 
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Due to the fluctuation of the heat noise, the equation of the elastica in heat bath obeys the 
MKdV equation (3-19-b) rather than the static MKdV equation (2-12). Let the set of solutions 
referred as 

©qci = {C\C is a solution of (13-19-b)}. (3-22) 

Of course 6qci contains ©cis- In Sqcb various shapes appear as the quasi-classical solutions in 
heat bath. For example, there should exist a deformed circle as a solutions of (3-19-b). As another 
example, there are other topological solutions of the different sector of the fundamental group, 

^ J d0qci G (3-23) 

Moduli of Closed Quasi-Classical Elastica. 

Here I will go back to compute the partition function (3-1), which was formally defined. First 

the problem aries how (3-1) should be interpreted. According to the philosophy of the canonical 
ensemble, the partition function should be calculated by the sum allover distinguishable and 
possible curves satisfying the isometry condition with Boltzmann weight. In the calculus, different 
topological class of (3-23) will be summed over. 

However the partition function (3-1) naturally diverges because the energy function is invariant 
for the affine transformation i.e. for the translation and the global rotation. Fixing C, if I denote 
X as 

X{s) =Xg+ Xr{s), J dsXris) = 0, (3-24) 

where Xg is the centroid of the curve C, the measure of the partition function can be rewritten 

as 

J DXe-f^^ = T.J ^^aj DXre-("'(^'> . (3-25) 
c 

The J DXg is volume of the base space CP^ and diverges. Including the rotational symmetry and 
the inner translation symmetry s — > s — , I will redefine the partition function which is divided 
by the volume of the affine transformation of the base space and length of the elastica, 

/ DXe-l^^ 

^■^^8 vol(Aff(CPi))L- ^^"^^^ 

Then I concentrate the shape of the elastica. I must classify the shape of the elastica and sum 
over the possible shapes. In other words, I must investigate the moduli space of the quasi-classical 
elastica, 

^'i^' Aff(Cp'i') X 51 • ^^"^^^ 

First I will consider the moduli space of the MKdV equation. The moduli of the MKdV 

equation was investigated as the KP-hicrarchy using Sato-theory [40-42]. (By the Miura map, 
the MKdV hierarchy is transformed to the KdV hierarchy and the KdV hierarchy is a subset of 
the KP hierarchy [41].) The moduli of the MKdV equation is classified with the genus 5 G N 
of the hypcrelliptic Riemannian surface (hypcrclliptic curve) Rg, which is the finite gap energy 
manifold (Bloch band spectrum) of the wave functions in the inverse scattering method of the 
MKdV equation [40-46]. Hence the moduli of the closed elastica is also classified by the genus, 

9^qci = n<cl- (3-28) 
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I will call the genus of the MKdV equation with the boundary condition genus of the elastica g. In 
fact, the classical solutions of (2-12) 's and (3-20) correspond to the elasticas of genus zero and one 
because these energy manifolds appearing in its inverse scattering method exhibit a Riemannian 
sphere and an elliptic curve respectively [43]. (It should be noted that even in the quasi-classical 
equation, the solutions of the circle and the eight-figure are unique up to homothety ioic g = 
and 3 = 1 respectively.) Using the knowledge of the properties of the universal Grassmannian 
manifold (UGM), I will consider the moduli of the closed elastica. 

First I will consider the simplest case {g = 0), a circle k = 27rn/L, n > 1. In (3-9), the 
quasi-classical action of these circles, 

^qci[C„] = ^, k{C^) = ^ (3-29) 

Hence for large n, the Boltzmann weight exp(— /3£^qci[C„]) rapidly decreases. This situation 
preserves for the elasticas of higher genus. 

Next I will consider elastica of genus one and why the number of the genus one solutions 
of closed elastica is only one, i.e., eight-figure shape, up to scaling. The moduli of compact 
Riemannian surface of genus one (or elliptic curve) is conventionally expressed as (1, r); r G TIr^ ; 

Mr, = H+/PSL{2,Z), 

H+ := {m e C\Im{m) > 0}, 

PSL(2,Z) := ad-bc=l,a,b,c,deZ^ /{±1). (3-30) 

However there is a dilatation freedom {K,K' := Kt) and thus I will denote TIr^ := M>o x 
to include the freedom: K G M>o := {x G M|x > 0}. The Jacobi variety of an elliptic curve 
is given as Ji,m := C/{KZ © K'Z) for m := {K,K'). Since is a real analytic function 
oi s £ = M/LZ, its domain embedded in Ji^m must be real. Thus only one-dimensional 
parameterization of (pqci^ C Ji^m, is allowed, which is direct line in Ji^m and passes its origin, 
because Ji^m is complex one-dimensional manifold. Since the moduli was divided by PSL(2,Z), 
there are PSL(2, Z) choices how such is embedded Ji.m- I choose a function with period L/n, 
n G N := {n G Z|n > 1} 0qci(s + L/n) = ^qci(s). By the periodicity of ^qci(s), I will fix the 
K for each embedding of into Ji,m) then the moduli of the MKdV equation with period L is 
OJIki X N xPSL(2,Z)/M>o, which is equivalent with N x 

On the other hand, the closed condition (2-19) in CP^ restricts the moduli of the elastica. I 
will introduce a real analytic map, 

a>t..xNxPSL(2,Z)^ 

h{m) = X{L)-X{<S). (3-31) 

Both if+ and CP^ are complex one-dimensional spaces. The moduli of closed elastica with genus 
one is given as the inverse image of zero point of /i, 

<ci = /r'(0)- (3-32) 

Due to the analyticity of the map /i, is zero-dimensional manifold. Thus the kind of the 

shapes are countable and due to the uniformity, there is only one-solution for each n G N. 

In ordinary computations [4,10-14,16,32,33], by reparameterizing above as R/{L/n)Z, one 
starts with C/((L/n)Z © K'Z), n G N and G H+ without dividing PSL(2,Z) and searches for 
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Similarly properties of the moduli of the closed elastica with genus g (> 1) will be investigated. 
It is well-known that by the Sato theory, the characteristic of the KdV hierarchy in the KP hier- 
archy is to characterize its energy manifold in the inverse scattering method as hyperelliptic curve 
in general (compact) Riemannian surfaces [40-42]. The Miura map from the MKdV hierarchy to 
the KdV hierarchy are bijective. Thus I will deal only with the hyperelliptic curves in this article. 
First I will denote the moduli of the hyperelliptic curve of g{> 1) as SPt^^. Its element is conven- 
tionally expressed as {Ig,Tg), where Ig and Tg are g x g matrices; Tg = (ri, • • • ,Tg) = (rij) and 
Ig = (ei, • • • , e^) is the unit matrix. As I did in = 1 case, I will deal with K{Ig, Tg), K G ]R>o 
rather than (Ig,Tg) itself. It is known that the dimension of the moduli of the hyperelliptic curves, 
^Rg, is 2^ — 1. Then I will also introduce a real 25-dimensional lattice for a point of the moduli 
m eMn^ := R>o x [44-46], 



uijKej + rijKTj (m^, rrij G Z) 



(3-33) 



and the Jacobi variety Jg,m := C^/Fm- If I determine a point m of ^Rg, I can uniquely con- 
struct the Jacobi variety, Jg,m '■= C^/Fm. From the soliton theory, if the coordinates of Jg,m as 
real manifold are expressed by t^^p = {ti,t2,ts,t4, ■ ■ ■ ,t2g) [40-42], its subset with odd indices 
tg '■= (^1)^3; ■ • ■ !^2,9-i) can be identified with the part of t in (3-12). This identification can be 
guaranteed by the Krichever construction of the solution of the KP hierarchy [44] and Sato theory 
[40-42]. By the Krichever construction, it is known that each parameterization G txp is direct 
line passing the origin in the Jacobi variety. (3-17) and (3-19) are reduced to the linear differential 
equations in the Jacobi variety [44-46]. (Thus (3-17) can be recognized as the Jacobi equation of 
the Jacobi-ficld of the system (3-14).) 

Since the moduli of the hyperelliptic curves has been also divided by a discrete group Sp{g, Z) 
[44] like PSL(2,Z) oi g = 1 case, for fixing m G ^Rg, there are N x Sp(cj',Z) ways to embed S^, 
as a period of (j), into Jg,m- The number of the ways are equivalent to cardinal of N x Sp{g, Z) as 
a set. As I impose its periodicity of L/n (n € N) on S^, the dilation parameter K is determined. 
Let such moduli space denoted as 

^Rg,s^ ■■= X N X Spig, Z)/M>o. (3-34) 

By choosing a point of the moduli space Tlji^^s^j g x g lattice Fj„,si is uniquely determined 
and Jacobi variety is given as, 

^Rg,s^ ^ {r^,si }, Jg,m,S^ := CVF„,si . (3-35) 

I will fix a point of the modulus m G 971^^ ^51 for a while. From the properties of the MKdV 
hierarchy, is a real analytic function of tg. As I can expand it around a point tg using the 
properties of its real analyticity, 



M^'9)= E (^no,■■■,nAt[-tl^it'^-ts^■■■it'2g_^-t2g-lr\ (3-36) 
no,--- ,ng 



tg must be a system of real parameters in Jg^ra,s^ ■ By analytic continuation, tg can be locally 
complexfied. On the other hand, the Jacobi variety has a canonical complex structure, 

J '■ Jg,m,S'^ Jg,m,S^, = —1, (3-37) 

which consists with its affine (vector) structure. By the structure J, there are set of the real 
^(-dimensional submanifolds {S^ ^ 51} which includes the orbit of s(g S^) as its one-dimensional 
submanifold. Due to the analyticity of over Jg^m,s^ > these complexfication of (3-36) cannot 
contradict with this complex structure J^. Then can be regarded as an element of ^g^m,s^- I 
will refer such embedding. 



12 



and the set of ao is expressed as Emo(Sc,,m,s'i ) Jg,m,s^)- Since (^3, • • • , i2g-i) need not be periodic, 
this embedding is measurable and Emo(Sg ^ 51, Jg^m.sO can be regarded as a subset of the 
Grassmannian manifold, Emo(Sg,m,si5 Jg,m.,s^) C GL(R, 2g - 1) /GL(E., g)GL{R, g — 1). 

Then I can construct the fiber structure for Emo(Sg rn,si ) <^g,m), because for a way to such 
embedding, there is the trajectory space ^g,m,s^ /S^ as fiber space, where {t^,--- ,i2g-i) G 
^g,m,si/'S'^- I refer the fiber bundle as 

TTtraj : ^ Emo(Sg,m,S'i , Jg,m,S^)- (3-39) 

This fiber space also depends upon the point of the moduli space of the Riemannian surfaces 
DTtijg^gi. Hence the moduli of the periodic solutions of <?I)qci(s, tg), which is written as 9?tperiod' 
has also fiber structure, 

7rperiod:2nSiod^2Jlfl^,s- (3-40) 

For each point m G 9JtR_^ 51, the fiber bundle (3-39) stands up as a fiber of 9?tperiod- 

By the closed condition, I must restrict the moduli space. Here I will consider a real analytic 
map like (3-31), 

fg '■ ^pel-iod CP"'', 

/,(m) = X(L)-X(0). (3-41) 
Consequently I obtain the moduli of the closed elastica, which is expressed as 

<ci = 4-'(0)c9Jt(l„d- (3-42) 

Since image of fg is real two-dimensional manifold, for g > 2, dim(9Jlq^|) > 1 and 9?lqci 
measurable. 

For simplicity, I will introduce the notations: 



n 



for (n, m) G DJl^f =o> := tt'^ (n, m) (3-43) 

where X^^^ is countable part of SJi^^^o ^t^=o,n' restriction of DJI^^^^q ^ point n G X^^^ , 
is measurable part. Here Sm,n has coordinate (fs, • • • ,t2g-i) =■ ^m- Hence there is a map from 
the moduli to the shape of the elastica. 

h : 3 (n, m, t„) ^ ci"^ (t^"^ ) C CP^ (3-44) 

For a perturbative deformation like distortion to an ellipse from a circle, the energy manifold 
in the inverse scattering method has infinite gaps (or genera) in general. Thus such deformation 
is expressed in the moduli of ^ 00 and due to the integrability of the MKdV equation, the 
deformation can be predicted like harmonic oscillator around a stable point. This picture is 
supported by the linearized method of the nonlinear equation and is also built in the above 
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Partition Function. 

As I finish to classify the solution space formally, I will consider the fluctuation of the elastica 

(n] 

again. It should be noted that there is above limit of the sequence U2 corresponding to the genus 
of the elastica. If I encounter for n = AT 

nui^^ = Xui^\ AgM, (4-45) 

like (3-20), then 5t2(N+m)+i oc St2N+i for m > because of (3-12) and 

dt.,.,r.,,A'"^"'^ = KA""^"^^ = fe.ci^^-nf ) = (3-47) 

Accordingly there needs no other fluctuation parameter n > N because these fluctuation vectors 
are linearly dependent. The sequence of (3-9) should be truncated according to the philosophy 
of the canonical ensemble. Thus I will denote such minimal integer, which is a function of the 
solution, as 

indo : C N{C) G Z. (3-48) 

However from the soliton theory [40-42] and above argument, for C G ^qcp I conclude that 
indo(C) = g. Avoiding meaningless divergence, I will replace the infinite series in (3-12) with 
finite sum from 1 to g depending upon the shape of elastica. 

Since the direction of Sti is along the tangential direction of the elastica Cqd, its effect has 
been treated as the integral of 5ti oc s in (3-26). On the other hand, 6t2n+i {n > 1) includes 
the normal direction fluctuation and I must integrate the Boltzmann weight over 5t2n+i space 
depending upon the genus of the elastica. Linear independence of these bases are guaranteed by 
above truncation. 

Then for a curve C G 971^^] , the heat fluctuation of higher order is expressed as 

I n „ n 

6^^^E[C,dt2m+l]= E -7^]l'^^t2rn.+^)J dsll{dt,^^^,)k^. (3-49) 

Here mi = part vanishes due to the periodicity 

I dsds (ll {dt,^^^,) k'^ =0. (3-50) 

On the other hand, if the set {rrii} does not contain rrii = component, the integral commutes 
with these derivatives, 

» n n „ 

y n = n j (3-51) 

i i 

Since J k'^ds is invariant for the time t2n+i {n > 0) development from the soliton theory, (3-51) 
vanishes. Hence I obtain that all higher order fluctuations vanish, 

d^'^^ElC, 5t2m+i] = 0, for n > 0. (3-52) 

In other words, the effect of heat fluctuation is given only through the energy functional of the 
quasi-classical motions. Since for a quasi-motion of genus g, the curvature is precisely given as 

kis.U.U-- - .U„_A = k„„As.U + -^^tn.U + -l^dto.--- .U„_,-^SU„_,). ('3-53) 
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I must integrate the Boltzmann weight cxp(— /3 J dsk'^) over all 6fs except 5ti. Using the trans- 
lation symmetry and freedom of the integration variable, I regard that t2n+i = 5i2n+i/v^- 

Consequently I obtain an explicit form of the regularized partition function (3-26), which is 
expressed by 

ZreM=Yl E (exp(-/?£^qcl[C])) 

1 

= Y.Y1 (exp(-/3^qci[C"])) 
g=o nex° 



+ 



EE / dm / (fl dhn-i] exp [-pE^,,{ct\t[::^)]) . 



(3-54) 



This is the exact form of the partition function (3-1) of the non-stretched elastica without diver- 
gence. In the second term, there appear the integration of the type of / dxe"^-^*-^^ Thus it is 
expected that the prefactor of the second term begins with the negative power of /3. For large 
/3, the second term is less than the first term. Hence for the zero temperature limit /3 — oo, the 
second term disappears and only the contribution of the genus and 1 survives. Noting that 
the moduli of the quasi-classical elastica with < 1 is equivalent with the that of the classical, I 
obtain 



lim ^reg[/3] = max cxp(-/3Ecis[C]) 

= exp(-/3 min Eas{C\). (3-55) 

Depending upon the boundary condition, the classical solutions appear as minimal points of the 
partition function Z^^^^\. Hence this partition function (3-54) does not contradict with the 
discovery of Daniel Bernoulli [1] . 

§4. Discussion 

It is worth while noting that due to the isometric condition, I have derived the MKdV hierarchy. 

In the clastic body, the lagrangian coordinate system should be employed rather than the Eulerian 
coordinate system when I will use the terms of the fluid mechanics. In the elastic body theory 
by marking some points on an elastic body and by estimating variation of distance among the 
marking points which is measured using the induced metric, the force will be locally evaluated 
as linear response for its certain deformation. The marking points corresponds to the Lagrangian 
test particle in the language of fluid mechanics. On the other hand, as I have used the metric 
induced from the base space and calculated the deformation, my calculation corresponds to the 
Eulerian one. Here it should be noted that if one uses the induced metric or Eulerian picture, 
any stretching (physical) curve can be regarded as a non-stretching (mathematical) curve; it is 
a trivial trick between the lagrangian picture and Eulerian picture and such recognition has few 
physical meanings. If stretching has physical meaning like an elastic body, Eulerian picture does 
not exhibit dynamical situation and unless stretch plays important role like boundary curve of 
binary fluid, dealing with stretch has less physical meaning. Accordingly the isometric condition I 
employed plays the central role in this scheme. In other words, in above computation, the reason 
why I could physically use the Eulerian picture even in the elastic body problem is owing to this 



15 



It should be also noted that even though there appears non-linear differential equation in 
this scheme, I have used the energy functional which is locally given in the framework of linear 
response of the force for the deformation [4,11]; if one uses the non-linear energy functional, he 
must evaluate it from basic elastic body theory because it is beyond the ordinary elastic body 
theory. It is remarked that due to the bilinearity of the energy functional, which is established in 
the framework of the ordinary elastic body theory, I could find the exact partition function (3-54) 
in this model. 

Furthermore the origin of the MKdV hierarchy in ref. [8,9,33] was artificial and was not phys- 
ically supported. If one physically sets up problem of time development of the elastica for real 
physical time, he concludes that its motion is not governed by the MKdV equation nor the MKdV 
hierarchy in general [10-14]. However in this article, I obtained the MKdV hierarchy from the 
physical requirement and a (mathematical) parameter time St2n+i appears variational direction 
as I pointed out in ref. [10]. In other words, by virtue of the novel investigation of the properties 
of isometric curve of Goldstein and Pctrich [8,9], I conclude that the virtual dynamics is realized 
as thermal fluctuation of an elastica in heat bath. Due to the isomctry condition, these equations 
become non-linear differential equations. In the linear differential equation such as the harmonic 
oscillator, the mode, which is determined by the global feature of the system, is represented by a 
vector of momentum space. As well as mode analysis of the linear system, these parameters t2n+i 
exhibit the global deformation of the elastica due to the thermal fluctuation and is expressed as 
a vector in the Jacobi variety. 

It is remarked that the obtained partition function (3-54) differs from that in ref. [24], which 
is obtained by summing the weight function over the conformation including non-isometry de- 
formation. Due to the isometry condition, non-linear terms appears in the quasi-classical curve 
equation while the partition function proposed by Saito et al. [24] is essentially linear. However 
for perturbative deformation, e.g. from the circle, the non-Iincar term might be negligible. Thus 
as long as the deformation is in perturbative, their partition function can be applicable even for 
a polymer which cannot be stretched even in thermal fluctuation. 

On the other hand, my partition function is justifiable even for large deformation. The partition 
function is summed over different topology which is related to the writhing number of the 
conformation. Hence there is a possibility of the topology change due to the thermal fluctuation. 
It is of interesting to calculate the possibility (or kernel function) from g = 1 conformation to 
g = 2 conformation. Even though the partition function (3-54) has not been concretely calculated, 
such computation, in principle, can be performed. 

Next I will comment on the physical meaning of 6t2n+i and the relation between the BRS 
transformation [31] and the Sato coordinate [40-42]. Since I have dealt with the SO (2) principal 
bundle over S^, the gauge group is expressed as 

C ]J S0(2), (4-1) 

where ]J means the disjoint union. (5 is infinite dimensional Lie group. It acts upon the shape 
of the curve, which corresponds to a section of the principal bundle, and deforms it. For a 
given shape of elastica, there is a unique group element which acts the elastica to become the 
shape with constant curvature, i.e., the simplest classical solution with g = 0. Thus the genus 
is well-defined, which is induced from the genus of curve (the quasi-classical section). There is 
a decomposition of as a family of subgroups <3g respect to the genus, whose action on the 
elastica preserves its genus. The representation of each group (Bg will be realized as (8g module in 
the set of corresponding Jacobi varieties. However in the soliton theory, instead of dealing with 
individual sets of Jacobi varieties of genus g, it is natural to consider the UGM if one wishes to 
formally treat a soliton equation. In fact there are singular elements in (3, which change the genus 
of the elastica; the transformation are known as the global gauge transformation. Corresponding 
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g and then naturally contains the singular elements due to the natural extension of t he group 
action. Thus is represented as a subset of GL(oo) in the UGM. The quasi-classical curve (a 
section of the S0(2) principal bundle) is embedded in the UGM. The infinitesimal deformation of 
the curve in the UGM can be expressed by a vector in the UGM. In other words, such deformation 
exhibits (mathematical) velocity of a trajectory in the UGM and can be represented as subset of 
the infinite dimensional general linear Lie algebra gl(oo), which is known as the afHne Lie algebra 
yl^^^ [41]; A^^'' is the Lie algebra associated with the Lie group (8. I will introduce the extrinsic 
differential operator in the UGM 

6:=Y,dt2n+idt.^+,, 6^ = 0. (4-2) 

n 

Then (3-15) and (3-17) are expressed as for A := fcqdds, 

SA = nui, (4-3) 

where 

"1 = "i"^dt2„+i, n := nk-J^d,. (4-4) 

n 

Noting the fact that u^"-* is hamiltonian density of the MKdV hierarchy, 

Sui K 0, (4-5) 

where « means equivalence after integration the both sides over s like (3-50) and (3-51). Since 
Ui obeys the Grasmannian algebra, Ui can be regarded as fermionic field over S^. Consequently 
(4-3) and (4-5) can be regarded as the BRS transformation of this system. Hence 5t2n+i in the 
path integration may be naturally understood in the framework of the Faddeev- Popov integration 
scheme [31]. In fact the square root of the Frcnct-Serret system (2-2) can be regarded as the Dirac 
operator [11,12], which is realized by confining the free Dirac field into a thin elastic rod. Using 
the Dirac field confined in the elastica, I constructed the MKdV hierarchy and r-function as the 
partition function of the Dirac field [47]. Thus I expect that the partition function (3-54) should 
also be expressed by the r-function of the MKdV hierarchy. 

Here I will mention a conjecture associated with the critical phenomenon of this elastica model. 
The critical point must be determined as a topological discontinuity of the moduli space of the 
quasi-classical elastica. At the point, physical quantity sometimes diverges and becomes meaning- 
less. I expect that the length of the elastica becomes less important around the critical point, e.g. 
the topological change, as a kernel function at an ordinary second order critical point becomes 
scale-invariant [48] . I will consider a dilatation of a quasi-classical elastica for the normal direction 
of elastica, 

Xc = X + ite'*% (4-6) 

where Xc and (pc are an affine vector and tangential phase of the quasi-classical elastica at the 
critical point and t is a real deformation parameter. The infinitesimal length of the elastica dsc 
becomes 

dsc = VdX^dX^ ^{1- kct)ds, := a.^c, (4-7) 
and then the length of elastica is 
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Noting that (pdL) — (pd^) = 27m, ra G Z, this deformation makes the length of the elastica 
change. However this deformation must be compatible with the isometric condition since I have 
been dealing with the non-stretching elastica. Both requires seem to contradict with but the 
critical point is an irregular point at which the contradicted objects coexist. From (3-2), it in 
(4-6) must be proportional to Uc ■= «i + iu2 with (3-18-b), and ui satisfies (3-19-b) in respect of 
the deformation parameter t. 

ict = Uc = - ^dskc, (4-9) 

where c is a real proprtional constant. Since this relation is the Miura map, from (3-19-b), Uc 
obeys the KdV equation 

ic = dfUc = GucdgUc + d^Uc- (4-10) 
By integrating it in s, (4-10) becomes 

ics = 3ul + \dluc. (4-11) 
For z ■.= \s w := Uc/2 and c = —2, (4-11) can be rewritten as 

d^^w = + z. (4-12) 



This is the Painleve equation of the first kind. Thus I will conjecture that at a critical point 
of the elastica, an expectation value obeys the Painleve equation of the first kind [49]. In our 
scheme, I can formally obtain the series of the ordinary differential equations related to the KdV 
hierarchy (from the Miura map) and the Painleve transcendents. Thus the muduli space must be 
closely related to the quantum two-dimesnional gravity, in which (4-12) and the KdV hierarchy 
naturally appear [50-52]. In fact, the loop soliton partially appears in the immersed surface in 
three dimensional space [38,53]. By the fermionic study [47,53-57], the immersed surface 
system is interpreted as a natural extension of the elastica system and also that of the Liouville 
surface system whose quatum version is known as the quantum two-dimesnional gravity [55-56]. 
Thus I plane to investigate the immersed surface system and reveal the relation between the 
elastica system and the quantum two-dimesnional gravity. (Here it should be noted that elastica 
problem is not directly related to the string problem in the string theory because the action of 
the elastica is biharmonic for X while that of string is harmonic [55]. In non-relativistic space, 
thickness is more important and biharmonic equation is very natural.) 

As I dealt with the kinetic properties of the large closed polymer and investigated the moduli 
of the MKdV equation in CP^, recently another statistical mechanics model of a large polymer 
was reported [58]. Partition function of non-contractible self-avoiding two-dimensional polymers 
in the topological torus, 

r = C/(LiZ + L2Z), (4-13) 

was studied associated with the MKdV equation. The partition function of such polymers was 
also solved by the MKdV equation [58] and is recognized as its r-function [59]. As I assumed the 
base space as CP^, 1 can also replace it with the topological torus. Then the ratio L : Li : L2 
becomes important as the boundary conditions but I can formally calculate it. Then obtained 
partition function may be also closely related to the r-function of the MKdV equation. It is a 
very interesting fact that the MKdV equation appears and plays central roles in both theories 
even though they are not directly concerned with as models. 

Finally I will mention farther possibilities and development of this theory. I investigated the 
elastica in two-dimensional space but can extend my theory to that in three-dimensional space 
if I can classify solutions of quasi-classical curves of elastica in three-dimensional space. Thus I 
need investigation on the moduli of quasi-classical elastica in three-dimensional space including 
topological properties like a knot invariance. Furthermore even though I formally classified the 
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of the closed elastica oi g > 1 now. Accordingly it is very important for this study to find explicit 
shapes of closed elastica of g > 1. If I could explicitly draw any shapes of closed quasi-classical 
elastica in plane of = 2 or 3 using hyperelliptic functions, I think, they have enormous effects 
on this problem. Since the summation of the partition function might be converges on the genus 
g, to determine the shapes even for small g means important steps on this problem. Moreover, it 
is also of interests to deal with an elastic rod which can be stretched as more physical model [14] 
and higher dimensional objects, such as an immersed surface [38,53,56,57]. 
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